“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1966 


An Automatic-Throttle Device for Airspeed 
Control During Flare to Field Landings. 


Elich, Robert Louis 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/30984 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


AN AUTOMATIC-THROTTLE DEVICE FOR AIRSPEED 
CONTROL DURING FLARE TO FIELD LANDINGS 


ROBERT L. ELICH 





L® 3m Aggy 
. oe ZL a 


WAJAL POSTGRADUATE SCHaOr 
MONTEREY, CALIF. 9394" 

















| DOWNGRADED 
APPROVED FOR PUBLIC RELEASE 


AN AUTOMATIC-THROTTLE DEVICE FOR AIRSPEED CONTROL 


DURING FLARE TO FIELD LANDINGS 
by 


Robert Louis plich 
Lieutenant Commander, United States Navy 
B.S., University of Colorado, 1957 


Submitted in partial fulfillment 
for the degree of 


MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 
from the 


UNITED STATES NAVAL POSTGRADUATE SCHOOL 
May 1966 


ABSTRACT 

This study is a continuation of the investigations in the use of automatic- 
throttle compensation systems. The previous investigations demonstrated 
the capability of such systems in alleviating the problems of aircraft speed 
instability in the carrier landing approach. The feasibility of adapting the 
system with sensed-input variables of angle-of-attack and normal accelera- 
tion to the condition of field landings with flare was investigated. An analy- 

oo 

sis of the ground effect prevalent in fiet@Tandings was performed and the 
significant influences on aircraft response were included in the investiga- 
tion. The desired airspeed control with zero rate of descent at touchdown 
was obtained by incorporating additional commands in the basic auto-throttle 
device during the flare maneuver. The analog computer was utilized for 
problem simulation. 

This investigation was performed by LCDR Robert L.iElich, USN, at the 


U. S. Naval Postgraduate School, Monterey, California. 
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Me Introduction. 

This study is a continuation of the investigations being conducted at 
the U. S. Naval Postgraduate School on automatic-throttle compensation 
systems. The previous studies conducted by Lieutenant Commander G. R. 
Bell [1], and Lieutenant R. E. Evans and Captain R. H. Schuppe [2] were 
presented in 1963 and 1964 respectively. 

As shown in Ref. 1 and Ref. 2 automatic-throttle compensation system 
utilizing various sensed-input variables are capable of alleviating the speed 
instability difficulties of modern supersonic aircraft during carrier landing 
approaches. One such system utilizes angle-of-attack change and change 
in acceleration normal to glide slope for commands to initiate the proper 
change in engine thrust. The purpose of this investigation was to determine 
the feasibility of adapting this system to the condition of normal field land- 
ing which is terminated with a flare maneuver to the touchdown. It is pro- 
posed that by incorporating additional commands to the basic auto-throttle 
system during the flare maneuver, a decrease in airspeed to a desired min- 
imum at touchdown can be accomplished simultaneously with zero rate-of = 
descent. 

The linearized, full equations of motion derived from small perturba- 
tion theory were used for this analysis. The coefficients in the equations 
were assumed constant except for those directly affected by the aircraft 
flying in close proximity to the ground. The influence of ground effect on 


the dynamic stability due to the varying coefficients was investigated and 
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was included in this study. 

The vehicle used in this and the previous studies is the Navy swept 
wing fighter, F-8. The Donner 3100 Analog Computer with associated equip- 
ment was utilized to simulate the aircraft motion and.maneuvering response 
with the proposed auto-throttle device. 

De Basic Automatic-throttle ( c<- n) System. 

The type of automatic-throttle controller presently used in the Navy 
F-8 is Specialties, Inc. Automatic-Power Compensator, (APC), which utilizes 
sensed-input variables of angle-of-attack, ©, and normal acceleration, n, 
for engine thrust command. A block diagram, typical of this system is shown 
in Figure 1. For the design application of this system, constant glide slope 
and airspeed are desired throughout the approach and the two sensed vari- 
ables of instantenous angle-of-attack and acceleration normal to the glide 
slope in "g" units are the engine command functions. The use of these 
variables as the command functions for airspeed control is evident from the 
relationship: 


a 
L =<"pn'w = yf S Oy, 
n'W = SP VS CL & 
ox ZPNV' W YLex 
Nn’ = Const. (V2) 
x 


From this equality, and increase of velocity results from a positive normal 
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acceleration increment and a velocity decrease results from a positive angle- 
of-attack increment. Thus, the normal acceleration in "g" units and angle- 
of-attack are inputs which are converted into the thrust command by the 
system computer utilizing chosen gain and time factors. The resulting change 
of thrust from the engine acts as a forcing function to the airframe modes of 
motion appropriately changing the aircraft velocity. The APC computer in- 
directly controls the aircraft velocity through the equations: 


7 A 


Kea Ax + Kj4n' 


KpAT 


KL.AX+ Ki An! 


since AV 


therefore, AV 
where the K values are the appropriate shaping, filtering and simulation 
networks. 

The transfer function for the thrust command as specified by Speciak.-_, 


ties,. Inc. is: 
Ra KY KY 
——e 


este 7g. leds 
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where, 
- o.2 sec (throttle servo time constant) 


- 0.5 sec (angle of attack filter time constant) 


1.2 sec (engine time lag) 


N a Q aN 


- 1.0 sec (normal acceleration filter time constant) 


K.. - 570 lb/deg. (angle-of-attack displacement gain) 


1S 


WN 
x 
l 


190 lb/deg-sec (angle-ofattack intergral gain) 


Kf 1,100 lb/deg (angle-of-attack gain) 


wN 


A 19, 400 lb/"g" (normal acceleration gain) 

As evidenced from the transfer function, the effect of increasing the 
instantaneous angle-of-attack results in an opening of the throttle and an 
increase of thrust, and a decrease in thrust is obtained for a decrease in 
angle-of-attack. An increment of normal acceleration is seen as a negative 
element in the transfer function with the effect of decreasing the thrust for 
a postive "g" perturbation and increasing the thrust for a negative "g" per- 
turbation. Ref. 1, 2, and 3 support the validity of this auto-throttle device 
for use in carrier landing approaches. 

3. Discussion of Auto-throttle System in Field Landings. 

In the field landing phase, it is desired that just prior to touchdown 
the vertical kinetic energy be dissipated and the horizontal kinetic energy 
be a minimum contingent on the minimum safe airspeed. Thus, structural 
damage and wear and tear on aircraft components will be minimized along 
with the advantage of using shorter runways. From the pilots point of view, 
these results are attained by application of elevator to flare the aircraft as 
necessary to decelerate to the desired touchdown airspeed with zero rate- 
of-descent. 

It is proposed that the &% -n auto-throttle system can be utilized to 
accomplish the desired airspeed control by incorporating a supplementary 


circuit in the designed APC. For the purpose of this investigation, it was 
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assumed that the F-8 aircraft was making an Instrument Landing System (ILS) 
approach with a four degree glide slope on the recommended approach speed 
of 1.2 Voz = 139 Kt. commencing a flare at 100 ft. above the ground, approx- 
imately eight seconds prior to touchdown, The desired airspeed for touch- 
down was 1.05 Vo; = 123 Kt. or 209 ft/sec. 

Since it is mandatory that the auto-throttle system be unchanged up 
to the flare maneuver, it was decided that the proposed supplementary cir- 
cuit of the APC be separately energized at this point. This circuit would 
supply the auto-throttle device with additional angle-of-attack and/or normal 
acceleration signals changing the error signals received. Thus, in effect a 
new base for the auto-throttle system is created. The error signal changes 
shall be referred to as "biasing" throughout the remainder of this study. 

The method of applying the angle-of-attack bias was determined from 
the equation describing the change in pitching moment coefficient from 
straight to rotational flight: 


AG. = aoa a mp 47) + mag 


Assuming that the change in moment is zero and neglecting the change in 


moment due to rate of pitch, q, the resulting equation is: 


oh ae elas Vm 47) 
Om 


Pg ny: 


OF : 
MX 


U 





AX = Const. (A7)) 


Thus, the angle-of-attack bias was experimentally utilized as some percent- 
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age of 4 7). 

Since it has been previously determined that the average pilot pulls 
a maximum of + 0.2 g above the 1 g steady state condition during the flare 
maneuver [4]. the normal acceleration bias was experimentally utilized as 
a percentage of 0.2 g. 

4, Equations of Motion. 

The basic three degree of freedom longitudinal equations of motion _. 
based on body axes used for this analysis are shown below. These equations 
were derived from linearized theory based on small disturbances from a ref- 
erence steady state flight condition. Full derivation and definitions are 
included as Appendix I. The equations are shown in the dimensional form 
necessary for analog simulation. 

X FORCE: 


, Xx Xx 
u- =u u- ax +g cos@o = aS] nop 
m m 
Z FORCE: 
ay 


mUo 





eee Z x 
eto 8 = OS 


MOMENT EQUATION: 


M Ma - ™M . “1. . mM M 
& woe OX = eee =x eee Se Al + Sell 
os 8p B 5° go) ee 








The values of the dimensional stability derivatives and aircraft parameters 
as obtained from Ref. 3 are included as Table l and Table II respectively. 


or Ground Effect. 


The principal effects of the ground on a nearby aircraft are increases 
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of the slope of the lift curves of the wing and tail and a decrease in the 
cdownwash at the tail. These changes result from the constraint of the 
vortex patterns at and behind the wing. The effects may be approximated 

by superimposing an image vortex system, as shown below, beneath the 
ground having induced velocities of necessary magnitude to satisfy the 
boundary condition of no flow across the ground plane! This method of 
approximation utilizes only the trailing vortices as it is normally considered 
that the bound vortex patterns cancel each other. The trailing image vortex 
will decrease the downwash at the tail and decrease the downwash in the 


tail and wing planes which increases the slopes of the lift curves. 


a a 
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For a two dimensional lifting surface, the Kutta-Joukowski theorem 


states that the lift is equal to the density times velocity times circulation 


or 


L = oV [ 1bftt. of span 


landrews, E. J., Notes on Simplified Ground Effect Analysis. U. S. 
Naval Postgraduate School, Monterey, California 


ro 


Thus the local lift coefficient is: 


Cm iT 
5 V4 Cc 


- aoa 
c V 
For an airfoil with uniform span loading with the associated trailing vortices 


a distance Kb apart, where b is the span, the induced velocity at the aircraft 


centerline is: 


——— 
4Tl Kb/2 


Since the downwash, € , is defined as =a then 
[* 


a7 V Kb/2 


Circ 
4 TF Kp Z 


Therefore, the change in downwash with respect to the change in the local 
lift coefficient is proportional to os . Resolving the induced velocities 
from the image vortex system into the vertical plane at the aircraft centerline 


results in: 





bee [" Kb/2 
ant | + 2ay2 ea + 7 








| CcV SH 
W aS 
47 7 als 
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d= Le) Saas 
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déi uw Kb/2 
C 
: AP)? + (2a)? 


Therefore, in proximity of the ground the change of downwash with respect 


to local lift coetficient is proportional to the following: 








cin | we Kb 
OC], EP? + (20)? 


By setting up a ratio between coh out of ground effect and g€_ in ground 
oC} SC) 


effect, the following equality is obtained: 


(Sc 

25 G _ l 
Kby 2 

va Le 


Aerodynamic texts state that for the lifting line theory of finite wings 
the trailing vortex filaments can be assumed to be 0.8 of the span apart. 


Substituting this value for K in the above equation yields: 


a€ ). 
dC; a 1 
— "> nari 
a€ ) 1 Sola 
Yor 
The ratio of QE in ground effect to 1 out of ground effect is only depen- 
oC, ac, 


dent upon wing span to height above ground, b/d. Therefore, a single curve 
utilizing the coordinates of the ratio of downwash versus b/d ratio can be 
constructed. This curve is shown in Figure 2. 
For a finite aspect ratio wing, the lift coefficient is given as 
Cr = 27 (& - HX) 


C 
where &X; is the induced angle-of-attack which is equal to =the 





for ellip- 


Za 


tical lifting surfaces. Substituting this value into the above equation and 
taking the derivative with respect to & gives: 


Of, pra 
Po no.4 et 


By utilizing the same procedure as inthe downwash effect, the followingratio of 


slope of the lift curves imandout of ground effect was obtained: 


(2x) 

aaKIG | A + 

ss) - l 

eae A a 
OX ‘ 2| 1 + wr 


Fora particular aircraft with a wing of known aspect ratio, a plot of the ratio 
of the slope of lift curves versus b/d can be constructed. For the F-8 with 

an aspect ratio of 3.39, sucha plot was obtained and is included as Figure 

3. 

Utilizing the curves of Figures 2 and 3, the coefficients in the equa- 
tions of motion containing stability derivatives dependent on downwash or 
slope of lift curves of wing or tail were varied. The stability derivatives 
effected were Xx, Zy, Mc, andMq. The definitions of these derivatives 
in Appendix I were rearranged as shown in Appendix Ilimorder to use the 
above derivation. Figures 4a, 4b, 4c, and 4d graphically present the value 
of these derivatives versus altitude above ground. The values of the coef- 
ficients containing these derivatives were then obtained and are presented 
in Figures 5a, 5b, 5c, and Sd. 

The proximity of the ground on an aircraft causes a lowering of the 


drag polar, reduced drag level for lift attained. This effect was considered 


ZZ 


offsetting the increased drag resulting from the increase of angle-of-attack 
during the flare maneuver. Therefore, the drag was assumed constant 
throughout this analysis. 
6. Analog Construction for Complete Problem Simulation. 

The equations of motion in dimensional form rearranged for analog 


simulation with necessary scaling factor notation included are: 





ayu = -.0602a,u + fe aq & - 31.88ag0 
+ .00145a,, ,AT 


cits = Z, — 7 
ax X = -.001l3a,u + —Lay & + a.0 


om rs 
- 0595 ay 7) 


x - Mé : M 
Be = 1 7oex 0 Fa u a ooo 4 ae ae 


- .33968 é - .455x 10°°a, AT - 2.50487) 
The scaling factors are included in Table III and the analog:mech anization 
of the equations is shown in Figure 6. 

The automatic-throttle transfer function given by Specialties, Inc. 
was converted into analog format as shown in Appendix III with the scaling 
factors included in Table III and then mechanized. The analog mechaniza- 
tion is shown in Figure 7. 

For the complete problem simulation, auxiliary equations for generat- 
ing actual altitude during the approach, glide slope and acceleration normal 


to the glide slope in "g" units were derived as shown in Appendix IV with 


scaling factors also ncluded in Table III. The equations 
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were mechanized as represented in Figure 8. These equations in analog 
format are: 


Altitude: 


ayH = 3 Ho - oltad) [ay Uo + ay W,+ aw /a x? 
Glide slope: 


35 o = 99 BKM 


Acceleration: 


==: > 


ll 17 eee “e 
Ain ea By Uot au) A, 9 - ag & ) 


a |p 


To perform the flare maneuver, an analog mechanization for elevator 
control was constructed by applying a steady voltage through a potentiometer 
to an integrator amplifier circuit. Thus, by adjusting the potentiomenter, 
application of step input or variable slope ramp input of elevator to the air- 
frame could be obtained. The elevator circuit was.further controlled by an 
offfon function switch built into the Donner computer. 

Te Experimental Methods and Results. 

7.1 Ground Effect Investigation. 

Due to an insufficient number of available amplifiers on the 
Donner Analog Computer for the simulation required in this study, it was im- 
possible to include all the ground effect variables. Therefore, it was necess- 
ary to evaluate the influence of each variable on the dynamic stability of 
the basic airframe inorder that a simplification could be obtained. Prior to 
experimentation, it was concluded that ground effect acting on the derivative 
M « would be the most influential and therefore would be retained. 


yw AA D. ayebgew, .1 ind 2 bide bed. +4 ! teh a ' - ans i (le a oF 
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A dynamic stability analysis on the airframe equations of motion with- 
out forcing functions of thrust and elevator was performed. Using cate 
condition of u(0) = 25 ft. /sec. for the phugoid determination and an initial 
condition of (0) = 5 degrees for the short period mode of motion determina- 
tion, the following analog runs were obtained: 

(1) Time histories of & and u perturbations out of ground effect in- 
cluded as Figure 9. 

(2) Time histories of & and u perturbations with all ground effect 
variables included for heights above ground corresponding to b/d ratios 
1.0, 2.0, 3.0, and 4.0. Results are included as Figure 10. 

(3) Time histories of & and u perturbations neglecting ground effect 
on derivatives Xqy and Z, for the same altitudes as in (2). The analog run 
outs are included in Figure ll. 

(4) Time histories of & and u perturbations neglecting ground effect 


on derivativesXy, , Z,, andM,« for the same altitudes as before. The 


w! 
analog representations are included in Figure 12. 

Analysis of Figure 9 determined that the phugoid period is 34 seconds 
and that of the short period mode is six. seconds. Times to damp to : ampli- 
tude are approximately 30 seconds and 2 seconds respectively. With the 
condition of all variables included, the maximum deviation from the no 
ground effect analysis in the phugoid and short period modes was 0.3 sec- 
ondss: There was no noticeable change in the times to damp , ampli- 
tude. Analysis of the analog runs when holding X,Y and Zw constant 


showed no noticeable change of period or time to damp to ; amplitude of 
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either mode of motion from those obtained with all derivatives variable. 
When the derivatives My, 4w, and Xy were held constant, the phugoid 
mode period was 0.2 seconds longer and the short period’s 0.1 seconds 
shorter then obtained when varying all the derivatives. The times to damp 
to - amplitude were basically unchanged. 

Since the relative influence of the variable derivatives Xx, Zy, and 
Ms, on the dynamic stability was inconsequential, the decision was made 
to conduct the investigation of this study with these three derivatives in- 
variant to ground effect. The derivative Ma which was retained as a 
variable dependent on altitude was constructed ona function generator for 
analog use. The analog mechanization is shown in Figure 13. 

7.2 Design Auto-throttle System Operation. 

The normal auto-throttle system with and without ground effect 
was investigated to insure proper analog mechanization prior to commencing 
experimentation for the proposed system during the flare maneuver. Analog 
runs with auto-throttle operating for initial conditions of horizontal gusts, 
u(0) = +5 knots, vertical gusts, % (0) = 5°, and pitch angle 6(0) = + 4°. 

Figure 14a shows the response obtained without ground effect for an 
initial horizontal gust. The velocity perturbation is oscillatory damped with 
an initial overshoot approximately 45% of u(0) and was without a second over- 
shoot. Figure 14b is the response for an initial horizontal gust with ground 
effect and agrees closely with that of no ground effect. The thrust applica- 


tion is oscillatory in both cases with a damped wave motion with a peak4 T= 


26 


FZ00°lb: 

Figures 15a and 15b present the response obtained for an initial verti- 
cal gust with and without ground effect respectively. The oscillatory veloc- 
ity perturbation is rapidly damped with an oscillatory throttle application 
with a peak AT = 6, 300 lb. 

The analog results for the initial pitch angle condition without ground 
effect is presented in Figure 16. The velocity perturbation is postively 
damped with a slight oscillatory tendency which is always negative. The 
thrust application is smoother than in the other initial conditions and has a 
peak 4 T = 600 lb. 

The results obtained for the design application investigated correspond 
to the results obtained in Ref. 1 and are as expected. The analog mechaniz- 
ation was therefore considered correct. The numerical forms of all the equa- 
tions are as in Appendices I - IV and the analog potentiometer settings is 
included as Appendix V. 

7.3 Elevator Angle for Flare Maneuver. 

The elevator deflection required to perform the flare maneuver was 
approximated from the equation given in Ref. 5 relating the elevator angle 
per "g" for a constant acceleration maneuver. This equation is: 
is a Sinn “ly i 7 © amare 

inh CL 


eo CL 
Substituting in the required values for the non-dimensional derivatives and 





C 
ra tt) 
an average normal acceleration of 1.19, the resulting elevator deflection 


2/7 


was approximately - 2.7 degrees. 

Since a pilot will normally pull a maximum of + 0.2 g above the steady 
state 1 g during flare, the rate of elevator application was obtained by ad- 
justing the potentiometer in the elevator circuit until a peak "g" perturbation 
of 0.2 g was obtained. 

Thus, elevator was applied at the required rate until a total of - 2.79 
was obtained or touchdown occurred. 

7.4 Determination of Auto-Throttle Supplementary Circuit. 

In order to obtain an insight as to what commands the supple- 
mentary circuit shall supply the design auto-throttle system analog runs of 
a flare maneuver with constant thrust and with the design system operating 
were performed. The constant thrust response is shown in Figure 17. The 
velocity is seen to decrease 18 ft. /sec. during the flare maneuver. The 
response obtained while the auto-throttle was operating is shown in Figure 
18. The velocity decreases initially then attempts to return to the beginning 
steady state velocity as a result of an increase of thrust. 

The above results indicated that the additional commands prescribed 
for the auto-throttle supplementary circuit must effectively cancel the com- 
mands that caused the thrust to increase. 

The analog mechanization of the two biasing commands, as discussed 
in Section 3, is shown in Figure 19. To simulate real system operation, 
the same gains were used in the angle-of-attack bias, car and normal accel- 


eration bias, gp, as in the original APC system; however, an additional 
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potentiometer was inserted in each circuit inorder to control the amount of 
&Xp Or Jp Sent to the auto-throttle block. The biasing command could then 
be comprised of both Xp, and gp or either one of them. 

Analog experimentation for determining the bias commands was per- 
formed as follows: 

(1) Biasing only the normal acceleration loop through a range of 
Omz g to 0.01 g. 

(2) Biasing only the <& loop through the range of 0.01 AT) tow, A7). 

(3) Biasing both the & and “g" loops through their respective ranges. 

The responses for the normal acceleration biasing were initial decreas- 
es in. airspeed to a minimum value followed by an increase. The minimum 
velocity was never less than 224 ft. /sec. regardless of the amount of Sp: 
as an increase of thrust was always obtained. The aircraft's motion was 
a pull up with a variable minimum altitude dependent on amount of gp. A 
representative record of the typical response obtained is included as Figure 
20 which is the run with gp = - 0.02 g. 

The typical response obtained when biasing only the angle-of-attack 
loop was a smooth decrease in velocity to a minimum followed by a slight 
increase to a new steady state value. The angle-of-attack increased rapidly 
to a peak angle then decreased to a damped perturbation about some median 
angle. The thrust was decreased initially by a maximum ofthree-to.: -2 %95 
four thousand pounds depending on the amount of bias and then began increas- 


ing as the flare continued, As & , was increased, a greater decrease in 
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velocity and higher peak angle-of-attack were attained. With the described 
elevator application, an \Ap of 0. 80277 resulted in a smooth velocity de- 
crease of 25 ft. /sec. to the desired 209 ft. /sec. at the touchdown point. 
Touchdown occurred seven seconds after flare was commenced but the ver- 
tical kinetic energy was not completely dissipated as a slight vertical veloc- 
ity component was evidenced in the altitude read out. The peak angle-of- 
attack = 3. 8° occurred at an altitude of 20 ft. and the maximum thrust de- 
crease was 3000 lb. This presentation is included as Figure 21. 

The results obtained during the analog runs while biasing both the 
angle-of-attack and rene! acceleration loops indicated that the airspeed 
response could be more advantageously controlled. Also, the magnitude 
changes in thrust were less than those obtained in the &p condition. The 
peak angles-of-attack were also decreased along with their rate of increase. 
The airspeed increased slightly as the flare was commenced then smoothly 
decreased to a minimum value. This minimum velocity was followed bya 
damped perturbation about a velocity approximately four ft./sec. greater 
than the minimum. With the elevator application as previously stated, an 
optimum condition was obtained with &p=. 80577 and gp=.15g. This 
bias condition resulted in the smoothest decrease in velocity to a minimum 
209 ft. /sec. that occurred exactly at the touchdown point. The flare was a 
smooth manuever with essentially all the vertical kinetic energy being dis- 
sipated. The peak angle-of-attack = 3.3° occurred 6.5 sec. after flare 
commenced and 3.2 sec. prior to touchdown. The thrust increased initially 


to a peak of + 400 lb. then decreased to a minimum - 2200 lb. and finally 


increased to a positive 1000 lb. above the beginning steady state value. 
8. Summary of Analog Results and Recommendations. 

The analog experimentation disclosed that a desired reduction of 
velocity could be obtained during the flare maneuver. Two different supple- 
mentary Circuits to the design auto-throttle produced a required touchdown 
velocity with a simplified elevator application. These circuits were com- 
prised of (a) biasing only the angle-of-attack loop by a percentage of the 
ramp elevator application, and (b) biasing both the angle-of-attack and 
normal acceleration loops with a percentage of the elevator and a constant 
value of "g" respectively. 

The Cy¥p produced a smoother throttle response; however, it hada 
greater magnitude change of thrust. This bias also caused a more rapid 
™ change and a greater peak c& while not dissipating all the vertical kinetic 
energy. The energy could be dissipated by changing the elevator deflection 
but the desired airspeed response could not be attained. The circuit com- 
prised of biasing both loops had a less smooth thrust response with a very 
slight velocity increase initially. The flare maneuver was smoother with 
essentially all the vertical kinetic energy dissipated. The rate of & in- 
crease was improved and a less peak & was obtained. 

In general, the supplementary circuit entailing biasing both auto- 
throttle input loops produces better aircraft response even though the thrust 
application was notas smooth. Further investigation to obtain an optimum 


supplementary circuit is required when all the ground effect variables are 
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included along with the changes in the lift and drag coefficients resulting 
from the increased angle-of-attack during the flare.. An analysis in system 
operation with gusts during the flare maneuver and the associated elevator 
application is also required. 

It is recommended that the required further analyses be performed on 
@ hybrid simulator consisting of analog and digital computers. These addi- 
tional analyses could not be performed in this study as the Donner 3100 
Analog Computer was completely saturated. 

a Conclusions and Acknowledgements. 

1. The use of the auto-throttle device with sensed inputs of angle- 
of-attack and acceleration normal to the glide slope for field landings with 
a flare maneuver is determined feasible. 

2. The throttle could be controlled to produce a reduction of airspeed 
during the flare to a desired minimum at touchdown together with zero rate 
of descent. 

3. The desired results were acpomplished by incorporating a sup- 
plementary circuit, energized during the flare maneuver, in the design auto- 
throttle system. This circuit would be comprised of an angle-of-attack rate 
command to a maximum and a constant signal command of normal acceleration. 

4. As might be expected, system functioning appeared to be critical- 
ly dependent upon the rate of application of elevator during the flare. 

The writer wishes to express his gratitude to Professor E. J. Andrews, 
Department of Aeronautics, U. S. Naval Postgraduate School, for his tech- 


nical assistance and encouragement in the preparation of this study. 
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TABLE 1 


LONGITUDINAL STABILITY DERIVATIVE VALUES 


= = ED; 
Xy = 41.14 GEG 
lb 

X — ae ‘ 2 

Ww a. ft. /sec. 
Xx = - 2272.14 Nels 

lb 

> 4 = ? el ee, 

AT ee ft. /sec. 

=a |) = * 
Ly = 181.69 ft. /sec. 
lb 

ZL = = : eee 

Ww oo ft. /sec. 
47 = - 9522.04 lb. 
M, = Rs ae es) lb. -sec. 
Ms = - 3940.41 ift.-lb.-sec. 
M & = - 109,377.45 = ft.-lb. 
Mg = - 32,602.94 #ft.-lb.-sec. 
Mar 7 - 0.437 lb. -sec. 
M7 = 240, 342.91 ft.-lb. 
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TABLE II 


F-8 AIRCRAFT PARAMETERS 


Uy = 234 ft. /sec. Cp, = 0.0 

c.g. = 24% MAC CLE = 2.998 

W = 22,000 lb. CL y = 0.3897 

A, = 8.1° Cm, = 0.738 X 1074 

eee Bae Coy = 70-38 

a = s7cutte. Cms = ~ 0.55 

Sp = 93.4 ft?. Cheyy = 1.622 

b = 35.666 ft?. h = 0.24 

leas 14- 0Rdutt. hn = Om292 

Gre lle7 pete Ori pe, = 2.387 

Cy = 0.87 QE = 0.4772 
om | 

Cp = 0.197 A = 3.39 


Sp. = QO. 963 
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TABLE Ill 


ANALOG SCALING FACTORS 


Maximum Voltage 100 volts 


z 


_— * 


2 


@ 


a 


100 v. = 50 ft./sec. ; a, = 0.5 ft./sec.-v. © 


100 v. = 50 ft./sec*. ; ay =0.5 ft. /sec®..-v. 
100 v. = 10° =.1745 rad. ; ay =0.1°/v. =.001745 rad/v. 
Jy _ rad. , » £ de 
100 v. = 10°/sec. = .1745 sec, | aa =.l mec roy. 
100 v. = 20°= .349 rad. ; ag =.2°/v. = .00349 rad/v. 
= 20° 2 rad. .  a- = .20/sec-y. = rad, 
100 v. = 20"/sec. = .349 Se i ag =. 2°/sec-v. = . 00349 See 


100 v. = 20°/sec? = 349 tad» ; age. 2°/sec2-v. = ,00349 Tad. vy. 
100 v. = 20° = .349 rad. ; ay =.2°/v. = .00349 rad/v. 
1l00v. =29 7 4, =e eZug/v. 

100 v. = 300 ft./sec. ; ay, = 3.0 ft. /sec. -v. 

100 v. = 8000 lb. ; ayy = 80 lb. /v. 

100 v. = 28.659 = .5 rad. ; an = 28.65°/v. = .005 rad. /v. 


a 
100 v. = 2.0 /secé, : et = .02 /sec2-v. 


100 v. = 300 ft. ; ag = 3.0 ft. /v. 


36 


FIGURE 1 


BLOCK DIAGRAM OF AUTO-THROTTLE SYSTEM 
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APPENDIX I 
EQUATIONS OF MOTION 
The basic three degrees of freedom aircraft longitudinal equations 
of motion as derived from small perturbation theory. based on body axes are: 
X Force: 
[m$z- X,]u - wx, + [mgcos® - Xy$z]@ = X(t) 
Z Force: 
-uZ,, + [™$z - Zay|W ~ | (mU -} Zq) Sz |e + mg sin®= Z(t) 
Moment Equation: 
- uM, - WMw - wMw #5 QB - aM, = M(t) 
Rearranging, substituting in forcing functions of thrust and elevator, and 


neglecting minor terms Xq, Xx , Z, sin@, and Zar the equations become: 





q! 
X Force: 
u - Yu - Ca + gcos@@ - “AT ay 
where xX = a. 
Z Force: 
- muy Zu t & - K-38 - 7 af 
Moment Equation: 
-—Yu — MA x + 9 — —1§ = At + — 
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Figure defining signs: 


rey 
CAL 


Definitions of Stability Derivatives: 


kw = FU SCL - Cp 


Xjt = cos( X%,+ {) 

2, = “PU SwiSp + Cy ) 
ae = - PU Sw (Cle, + Cp) 
2m = -4 Ug SwOL 


= ] 
My = OU SwrlC, + 5U oS m,) 


— lg 12 dé 
Mx = Uo(- o»-qtly CLT dat? 
Max = UgMy = Ul U St Cm.) 


ie at = 
Mq = a CPU Swe Cmq) 
M, 7 me - 0. 437 


2 
M7 = + PU ds t Cry ) 
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Substituting the values given in Table I: 
X Force: 
u + .0602u + 3.325% + 31.886. = .00145 AT 
Z Force: 
OO1ISU + K + .4271K - @ = -.059577 
Moment Equation: 
-1.89 x 10°-4u + .041& + 1.1390 + 6 + .3396 6 


oa 0s AT = 2.504 7 
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APPENDIX II 


STABILITY DERIVATIVES FOR GROUND EFFECT 


ie Z 2C 
Kee ZO Ye Sw Gr (7 ~ 2 ous | 


(C.,. re 


ih 2 ses DN RAS 
2% Si Ga [ ~ eA oe (Cu. ) 


O€ 
= a, VY, CL ( a 


= Curie le 
= Cc. ( ~ Ps ae 
Cormlass yes:) Bele ($& \, 


Y —_ 
Mar (Gey [MARES (38) 





O€ 
YN bee OX 
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APPENDIX Il 


ANALOG FORM OF AUTO-THROTTLE EQUATIONS 





AT l i Kee a 4 Ky | Cg | Kn n 

(1+ 7,8) (1 +7.) lt 7 & leis 1+7 S* 

_ ] ] 570 190 1100 )}K, _ 19,400, 
(1+ 1.2S) (1 + . 28) 1+.5S S 1+8 ies 


Equation without time elements 

AT = 570(57.3)&% + 190(57.3)/ + (1100)(57.3)% - 19, 400n 
anToaT = 32,66la, K + 10,887 an, /O + 63,000 a, - 19, 400 ann 
AT = .71& + 237K + 1.371% - 4,857 

where: 


Ke = 5701b./fdeg. = 32,661 lb. /rad. 


Ky = 190 1lb./deg-sec. = 10, 887 lb. /rad-sec. 
K, = 11,000 1lb./deg. = 63,000 lb. /rad. 
K, = 19,400 lb./g 
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APPENDIX IV 


ANALOG DERIVATION OF AUXILIARY EQUATIONS 


Acceleration: 
T= bh (normal to flt. path) = Vrotal oO 
d 
— — - 
(U” # an) AR (0 ) 


= (Uo + u) (0 - &) 


ann = (ay Uo + ayt) (48 - ago) 

02H = (3.0 Uy + .5u) (.00349 6 - .001745 &) 
San ti a + .1670) (524 8 - .262&) 
g 100 ° 

Glide Slope: 
From Appendix I; O=eo -x 
ay x = a8 - ag & 
00349 § = .00349 6 - .001745@& 


Qa kk 


x 
\ 


Total Altitude: 


H, - Initial altitude = 300 ft. 
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apH 


3,08 


ie U,(t) sin & + h 
= Ho - % (rad) [ Uy + ui 
= Hop- Vem Um & 
= Ho - .07/ Wewete (epi 2 
an Ho - ree ay Uo + fay Uo at au) fa, ¥ 


3.0H, - .07 / 3.0U, + (3.00, + .su) f .00349 4 


. _—— 
Hi = 07 f Uy + (U.. tu: . 167 a)(. 0087) f 
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APPENDIX V 


ANALOG POTENTIOMETER SETTINGS 


Potentiometer Coefficient of (Equat. ) Setting 
] u (u) . 0602 
2 8 (u) 2252 
3 X (u) .0116 
4 x (x) , 4271 
5 u (&) . 3238 
6 u (6) 0271 
7 & (9) 0205 
8 6 (6) . 3396 
9 XK (8) ce 

10 Vr (h) ne7 
ll & (h) a6 

12 X(K, - T) owe 
13 n(K, - 1) . 485 
14 x (Ky - T) p37 
15 x(Ky - 1) wa 
16 Ip variable 
17 uta) . 093 
Dil 7) (8) . 358 
22 1) (&) 1705 
23 Uy Wy) yi 
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Potentiometer 


24 


29 


26 


28 


29 


30 


31 


32 


33 


34 


35 


36 


37 


38 


Coefficient of (Equat.) 
u (Vp) 
Up (B/S Up) 
Hy (H) 
@ (An) 
& (An) 
g (An) 
Xp (K) - 1) 
Xp Ky - T) 
Ap Ky, - T) 
Funct. Gen. 
AT (u) 
A T (8) 
~B 


7) 
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Setting 
‘Loy 
~0275 
a0 
924 

» 262 

7a 1 
«lov 
237 

ay 

al 

» 250 

. 1043 
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